Recurrence properties of quantum observables in wave packet dynamics 
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We investigate the recurrence properties of the time series of quantum mechanical expectation 
values, in terms of two representative models for a single-mode radiation field interacting with a non- 
linear medium. From recurrence-time distributions, return maps and recurrence plots, we conclude 
that the dynamics of appropriate observables pertaining to the field can vary from quasiperiodicity 
to hyperbolicity, depending on the extent of the nonlinearity and of the departure from coherence of 
the initial state of the field. We establish that, in a simple bipartite model in which the field is effec- 
tively an open quantum system, a decaying exponential recurrence-time distribution, characteristic 
of a hyperbolic dynamical system, is associated with chaotic temporal evolution as characterized by 
a positive Liapunov exponent. 

PACS numbers: 05.45.Tp, 05.45.Mt, 42.50.Ar, 42.50.Dv, 42.50.Md 



I. INTRODUCTION 



Poincare recurrences of classical dynamical systems 
yield a great deal of information on their ergodicity 
properties Recurrence-time statistics complements 
and augments the information obtained from other well- 
known signatures and quantifiers of dynamical behavior 
such as return maps, Liapunov spectra, time-series anal- 
ysis, generalized dimensions, and so on. 

There exists a body of rigorous results in the the- 
ory of dynamical systems that pertains to universal 
properties of recurrences. Most notably, for Axiom-A 
systems!, i, Q and more generally for uniformly hy- 
perbolic systems [il, the recurrence time to a sufficiently 
small cell in phase space is exponentially distributed. 
Moreover, successive recurrence times are independently 
distributed, so that the sequence of these recurrence 
times has a Poisson limit law. These arc familiar fea- 
tures of stochastic systems (e. g., aperiodic Markov 
chains [6j, lZ, 8J), but the significant point is that they are 
also exhibited by deterministic systems with a sufficient 
degree of mixing. More detailed studies of recurrence- 
time statistics have been carried out in the case of low- 
dimensional chaotic systems, in particular, in the frame- 
work of one-dimensional maps. These investigations en- 
able us to distinguish clearly between different degrees of 
randomness in the dynamics, ranging from quasiperiodic- 
ity through intermittency to fully-developed chaos [9l.lTo|. 
Departures from the standard results for hyperbolic sys- 
tems have been analyzed, showing (for instance) that 
intermittency leads to power-law (as opposed to ex- 
ponential) distributions of recurrence times, with cor- 
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relations between successive recurrences [llj and non- 
Poisson limit laws [l^l • A variety of results is also known 
for recurrence-time distributions for Hamiltonian sys- 
tems and for measure-preserving maps that model as- 
pects of such systems. Interesting universal asymptotic 
properties, including power-law distributions of recur- 
rence times, can arise here owing to the highly non- 
uniform nature of the invariant sets in phase space and 
the stickiness associated with the remnants of invariant 

tori[il B d, d, [13, Hi, [M iOi • 

Besides recurrence-time distributions, recurrence plots 
comprise a related technique for analyzing dynamics 
ranging from periodicity to chaos [111, [22] . Several distinc- 
tive features of recurrence plots have been identified and 
established as indicators of specific kinds of dynamical 
behavior such as multiple periodicity, intermittency, and 
chaos. Recurrence quantification analysis [l^] has been 
developed in an attempt to deduce quantitative informa- 
tion from the heuristics of recurrence plots. 

In contrast to these detailed results for classical dy- 
namical systems, the ergodicity properties of expectation 
values of observables in nonlinear quantum systems, re- 
garded as dynamical variables, are much less comprehen- 
sively understood. A fundamental difficulty that imme- 
diately arises here is that the 'phase space' is effectively 
infinite-dimensional, because the complete information 
contained in a specific quantum state can be obtained, 
in principle, only if the mean values as well as all the 
higher moments and cross-correlators of all the quantum 
operators pertaining to the system are included in the 
set of dynamical variables. Moreover, significant roles 
are played by the precise nature of the initial state of 
the system, and by the interaction of the system of inter- 
est with the environment. Owing to these factors (and, 
of course, the non-commutativity of different operators), 
the dynamics of expectation values may be expected to 
be quite complex, even in relatively simple systems. It 
is therefore relevant to examine specific tractable models 



2 



of quantum dynamics in order to disentangle the effects 
of various factors, and to discern systematic trends. 

The time evolution of quantum mechanical wave pack- 
ets provides an appropriate framework for the pur poses 
described above. The phenomenon of revivals |24l| pro- 
vides a manifest analog of recurrences in coarse-grained 
dynamical systems. An initial state |'0(O)) has a revival 
time Trcv if the overlap function | (?/'(O)|'0(T'rov)) P re- 
turns to an e-neighborhood of its initial value of unity. 
A physical realization of this phenomenon is provided by 
the propagation of a radiation field in a nonlinear medium 
under suitable conditions. A wave packet of the radiation 
spreads and loses its original form almost immediately 
owing to the field-atom interaction, but displays revivals 
at integer multiples of Trev Correspondingly, all expec- 
tation values return to the neighborhoods of their initial 
values. Distinctive signatures of the revivals, and of the 
frational revivals in between successive revivals, can be 
identified by analyzing the (projections of the) 'phase 
trajectories' of the system in subspaces of appropriate 
expectation valuesf25|. Typically, the behavior of these 
trajectories is similar to that of the phase trajectories in 
quasiperiodic dynamical systems with many incommen- 
surate frequencies. 

However, departures from this scenario occur with in- 
creasing deviation of the initial state from perfect coher- 
ence, or increasing nonlinearity of the medium, or both: 
It is found [26] that regular revivals of wave packets no 
longer occur under these conditions, namely, high non- 
linearity and non-coherent initial states. A more detailed 
examination of this physically interesting, and obviously 
more generic, regime is called for. This is the primary ob- 
jective of the present paper. The objects we study are the 
time series formed by the expectation values of relevant 
operators. The aspects we focus on are the recurrence- 
time or first-return-time statistics of these time series, 
and their recurrence plots. 

In what follows, we use two simple models that capture 
the essence of the effects we seek to clarify and highlight. 
In the first of these, considered in Sec. |TT1 we use the 
fact ^27*1 that certain aspects of the dynamics of the radia- 
tion field propagating in a nonlinear medium can be mod- 
eled by an effective Hamiltonian involving the field oper- 
ators alone, without explicitly invoking the atomic opera- 
tors. The purpose here is to understand the effects of the 
lack of coherence of the initial state of the wave packet 
upon its subsequent evolution and upon recurrence-time 
distributions. The Hamiltonian concerned has a purely 
discrete spectrum, which implies, according to a rigorous 
result ^2^], that the dynamics of expectation values can- 
not be chaotic (the maximal Liapunov exponent cannot 
be positive). Nevertheless, wc find that, with increasing 
nonlinearity and/or lack of coherence in the initial state, 
the recurrence-time distribution of a generic observable 
changes from a form that is characteristic of quasiperiod- 
icity to an exponential distribution, which is customarily 
associated with hyperbolic (and hence chaotic) motion in 
classical dynamics. An exponential recurrence-time dis- 



tribution is therefore not restricted to chaotic systems. 

While this model serves to illustrate the point just 
made, it must however be recognized that the radiation 
field is really one of the two components of a bipartite, 
interacting, field-atom system. In order to investigate 
what happens in this situation, we consider, in Sec. Illli 
a bipartite Hamiltonian that involves both the field op- 
erators and the atomic ladder operators. Focusing on 
operators involving the field mode alone then amounts, 
effectively, to considering an open quantum system inter- 
acting with the (atomic) environment. Even if the initial 
state is a non-entangled direct product state of the two 
modes, entanglement sets in during temporal evolution. 
The controlling parameter in this case is the ratio of the 
respective strengths of the nonlinearity and the inter- 
mode coupling. For low values of this ratio, an entangled 
state periodically factorizes into the initial product state 
(apart from an overall phase), and revivals occur. For 
higher values of the ratio, such revivals disappear. De- 
pending on how far the initial state of the field is from 
perfect coherence, the dynamics of relevant expectation 
values emulates that of a hyperbolic dynamical system. 
Recurrence-time distributions, recurrence plots, as well 
as a calculation of the maximal Liapunov exponent based 
on phase-space reconstruction using the time series of an 
expectation value, all corroborate this conclusion, which 
therefore has implications for the dynamical behavior of 
open quantum systems. 

II. SINGLE-MODE NONLINEAR 
HAMILTONIAN 

As stated in the Introduction, we begin with a simple 
effective Hamiltonian for a single-mode electromagnetic 
field interacting with the atoms of a nonlinear medium, 

= ft(xa^'a' + x'a^V), (1) 

where the photon annihilation and creation operators 
satisfy [a, a^] = 1. The first term in H is the standard 
one modeling a Kerr medium, with a coupling strength x- 
The second term is an additional cubic nonlinearity with 
a strength x' ■ It introduces a second natural frequency 
into the system. Since x ' Ix is arbitrary and generically 
not a rational number, it modifies the relatively simple 
behavior in the absence of such a term. Both terms 
are of course diagonal in the number operator a^a, since 
a^'^a? = a'^a{a^a — 1) and a'^a?' — a{a^ a — l){a^a — 2). 
As a consequence, the mean photon number in any state 
of the system remains equal to its initial value as the 
state evolves in time. 

We first recapitulate very briefly how revivals occur in 
the absence of the cubic nonlinearity, i.e., when x' = 0. 
At the instants t — -K/ikx) (where fc is a positive inte- 
ger), the time-evolution operator U{t) = exp{—iHt/h), 
which is diagonal in the Fock basis, has matrix elements 
exp [— J7rn(n — l)/k], where n denotes the eigenvalues of 
a^a. Even though the exponent here is quadratic rather 
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than linear in n, this exponential has interesting periodic- 
ity properties [29|] as a function of n. As a consequence of 
these properties, U{'n /kx) can be expanded in a Fourier 
series with exp {—2'nij/k) as the basis functions, where 
the integer j runs over the range 1 < J < (A — !)• If the 
initial wave packet is the Gaussian wave packet corre- 
sponding to a minimum-uncertainty coherent state (CS) 
\a) (where a\a) = a\a) and a G C), then the state at 
time 7r/(fcx) can be written as a finite linear combina- 
tion of the states jae"^''*-'/'^) (for odd k) or the states 
|q, gi7r/fc g-27r4j7fc^ (for cven k). But these are again co- 
herent states, with parameters that are just phase-shifted 
versions of the original parameter a. It is then easy to see 
that, at time Trov = tt/x (corresponding to k = 1), the 
initial state revives for the first time. Periodic revivals 
occur at integer multiples of Tj-ov Fractional revivals are 
also evinced at specific times in between successive re- 
vivals. 

The time evolution is quite different when x ' 7^ 0. For 
generic initial wave packets and (irrational) values of the 
ratio X '/ x, exact revivals do not occur. Correspondingly, 
in the space of observables, periodic returns of observ- 
ables to their initial values is replaced by quasiperiodic- 
ity. Now, the dynamical behavior of expectation values 
in this model cannot become chaotic for any parame- 
ter values (as already mentioned), because the evolution 
is governed by a quantum Hamiltonian with a purely 
discrete spectrum [28|. However, as the mean photon 
number increases, or the initial state departs from per- 
fect coherence, the quasiperiodicity involves an increas- 
ing number of incommensurate frequencies. Concomi- 
tantly, the recurrence-time distributions and recurrence 
plots go over from the forms characteristic of quasiperi- 
odity to those pertaining to a dynamical system with a 
much higher degree of mixing, as we shall see below. 

Our aim is to examine the role played by the extent 
of coherence of the initial state on the recurrence prop- 
erties. An obvious dynamical variable for this purpose is 
the expectation value of the quadrature x — (a -I- a^)/ 
The standard CS \a) serves as the reference initial state 
of the field mode. The mean photon number in this 
state is of course |ap, which we shall denote by v. 
The non-coherent initial states we consider are photon- 
added coherent states fPACS) [soj . The normalized m- 
photon-added coherent state \a,m) has a quantifiable 
(and, in principle, tunable) departure from perfect co- 
herence and Poissonian photon statistics. It is defined as 
\a,m) — (flt)™ \a) /[m\ Lm{—i')Y^'^ , where m is a posi- 
tive integer and Lm is the Laguerre polynomial of order 
m. The mean photon number in this state is given by 
(to -I- \)[Lm+i{~v) / Lm{—v)] — 1. Moreover, the vari- 
ance of the photon number can be shown to increase 
sublinearly with the mean photon number, implying that 
the photon statistics is sub-Poissonian. With the exper- 
imental production and characterization of a 1-photon- 
added coherent state by quantum state tomography [3l|, 
photon-added coherent states are expected to play an in- 
creasingly significant role in future investigations. 
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FIG. 1: First return time distribution F\ of the expectation 
value (a;). Initial state |a), with v = \. 



The procedure we adopt is straighforward. We start 
with a given normalized initial state ['(/'(O)), let it evolve 
unitarily according to the Hamiltonian i7, and calculate 
{il){t)\x\^\;(€)) = {x{t)) in time steps St. This produces 
a time series. The interval of values of {x{t)) is coarse- 
grained into small cells, and a recurrence to a given cell 
is marked. This enables us to build up the recurrence- 
time distribution for the cell concerned. The time series 
is also used to draw a recurrence plot. We have checked 
that the results reported below hold good qualitatively 
for generic cell sizes and locations in the phase space 
considered. For definiteness, we present in Figs. [lHS]the 
results for a cell size 6x = 10~^ and a parameter ratio 
x'/x = 10~^. The data set comprises a long time series 
of about 10^ points, and the time step 6t — 10~^. Time 
in units of St is denoted by r. 

We now recall that, for uniform quasiperiodic motion 
with two incommensurate frequencies, the distribution 
of the time of recurrence to any cell (on the 2-torus) 
has a support comprising just three values [Tol. [s^. This 
result is derived from a 'gap theorem' for interval ex- 
change transformations (specifically, the rotation map 
on a circle), and follows from certain number-theoretic 
properties of continued fractions 33, 34]. For quasiperi- 
odic motion with more than two mutually incommensu- 
rate frequencies (i. e., on an n-torus with n > 2), the 
recurrence-time distribution remains a discrete distribu- 
tion with support at a finite number of values (although 
no generally applicable formula is known for the actual 
number of such values, for arbitrary n). Returning to the 
problem at hand, we find that when the initial state is a 
CS ja), with a mean photon number v that is of the order 
of unity, the recurrence-time (or first return time) distri- 
bution Fl of (x) is a discrete distribution with support at 
a small, finite number of points, as in Fig. [1] If, instead 
of the CS \a), we use the PACS \a, m) as the initial state, 
the points of support of the distribution Fi increase in 
number with increasing to. Figure [2] illustrates the case 
TO = 5. As TO takes on larger values, this trend rapidly 
leads to an essentially continuous, decaying exponential 



4 



0.1 r- 

0.08 - 

0.06 - 

Fl 

0.04 - 

0.02 - 



20 40 60 80 100 

X 

FIG. 2: Fl for the initial state |a,5) with z/ = 1. 
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FIG. 3: Fl for the initial state |a) with i/ = 100. 



distribution. The same trend occurs even for an initially 
coherent state \a), provided the mean photon number 
f ^ 1, as seen in Fig. [31 The onset of an exponential 
distribution is even more pronounced for an initial PACS, 
as may be seen in Fig. 21 

The recurrence plots obtained from the time series for 
(a;) are completely consistent with this scenario. Figure 
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FIG. 4: Fl for the nitial state \a, 5) with u = 100. 




FIG. 5: Recurrence plot of the time series {x). Initial state 
|q) with 1/ — 1. 




FIG. 6: Recurrence plot of the time series of (a;). Initial state 
la, 5) with u = 100. 

O depicts the recurrence plot in the case of an initial 
CS I a) with v = 1. The regular, patterned structure of 
the plot is characteristic [22| of quasiperiodicity with a 
relatively small number of incommensurate frequencies. 
In marked contrast. Fig. [HI shows the recurrence plot for 
an initial PACS |a, 5) with i' = 100. It is evident that, 
at the very least, a significant degree of mixing has set 
in, corroborating what we have already deduced on the 
basis of the recurrence-time distributions. 

We are thus led to the conclusion that the generic dy- 
namics in this model makes a gradual transition from 
quasiperiodicity to mixing as the mean photon num- 
ber increases, or if the initial state deviates significantly 
from perfect coherence, or both. It has already been 
mentioned that hyperbolically unstable dynamics is pre- 
cluded in this model. We have checked this out inde- 
pendently by augmenting the foregoing with an analy- 
sis of the time series obtained for (x) along the lines 
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customary [H, [H, [sj in the study of dynamical sys- 
tems, namely: phase space reconstruction, the estima- 
tion of the minimum embedding dimension of the effec- 
tive phase space, and the calculation of the maximal Li- 
apunov exponent Amax- The result, as expected, is that 
Amax = to within numerical accuracy. The occur- 
rence of an exponential recurrence-time distribution is 
not inconsistent with this result. As is familiar from re- 
sults for Markov processes, exponential recurrence-time 
distributions and statistically independent successive re- 
currences with a Poisson law can occur even in sta- 
ble stochastic d yna mics [ssf with a (mean) maximal Lia- 
punov exponent [39| that is zero or negative. The essential 
requirement is an adequate degree of mixing. 

III. TWO-MODE MODEL: OPEN QUANTUM 
SYSTEM 

Having seen how the introduction of the cubic non- 
linearity in the single-mode model of Sec. alters its 
dynamical behavior, we now turn to a more representa- 
tive model of the radiation field propagating in a non- 
linear medium. This model involves both the field and 
atom modes explicitly, and is given by the bipartite 
Hamiltonian 

H = h[uj a^a + Loob^'b + -fb'''^ 6^ + 3(0+6 + (2) 

As before, a and a+ are the annihilation and creation op- 
erators for the field. The medium is modeled by an an- 
harmonic oscillator with ladder operators b and b\ the 
parameter 7 serving as a measure of the nonlinearity. We 
focus on the relative roles of the nonlinearity and the cou- 
pling of the two modes (quantified by the coupling con- 
stant g) in determining the dynamics of observables. In 
the numerical results presented below, we have therefore 
set the unperturbed frequencies cu and ujq equal to unity. 
The relevant tunable parameter is then the ratio "//g. 
The idea is to examine the expectation value of an oper- 
ator pertaining to the field alone, which may be regarded 
as an open quantum system in interaction with the 'en- 
vironment' represented by the medium. Since an open 
subsystem of a Hamiltonian system is not conservative 
by itself, its evolution is effectively guided by a positive 
map, and any hyperbolicity in the dynamics could also 
lead to chaotic behavior of the observable concerned. 

The total number operator (a'^a + b^'b) commutes with 
the Hamiltonian ^ for all values of the parameters in H, 
but the photon number operator a+a does not do so for 
any g 7^ 0. Thus, while H can be cast in block-diagonal 
form in a direct-product basis of field and atom Fock 
states, the model is not trivial. The Fock basis is given 
by I") field I"') atom' where n and n' are the eigenvalues 
of a^a and 6+5 respectively. The basis states can be con- 
veniently taken to be |iV — n)^^)^ (g) \n)^^ = \N — n;n), 
where N labels the eigenvalues of (a^a -I- 6+6). But 
{N-n;n\H\N'-n';n') = if TV ^ N' . Hence, for 
a given value of N, H can be diagonalized in the space of 



states |iV — n ; n) , where n = 0, 1, ... A''. A natural choice 
for the observable representing the field (the 'open' sub- 
system) is the photon number operator a+a. 

It is evident that if g = 0, is just the sum of two 
decoupled parts. If 7 = 0, if is linear and hence is again 
trivially diagonalizable, and there is periodic exchange 
of energy between the field and the atomic oscillator. 
When 7 and g are both non-zero, the dynamics is more 
complicated, but the existence of the conserved opera- 
tor (a+a -|- 6+6) ensures that the system as a whole is 
well-behaved, with a discrete spectrum labeled by the 
quantum numbers N and s, where A^ — 0, 1, . . . and 
s = 0,1,..., A^. The time-dependence of (a+a) is thus 
a direct consequence of the coupling of the field mode 
to another degree of freedom; and the deviation of (a+a) 
from periodic temporal variation is an indicator of the ef- 
fects of the nonlinearity in the second degree of freedom. 
For different ranges of values of the parameter ratio 7/g, 
a diversity of temporal behavior is exhibited by (a+a). 
In the weakly nonlinear case (7/(7 ^ 1), collapses and 
revivals of (a+a) occur almost periodically in time for ini- 
tial field states which arc Fock states or coherent states. 
When 7/5 ~ 1, such collapses and revivals occur more 
irregularly if the field is initially in a coherent state, as 
compared to an initial Fock state. In the strongly non- 
linear regime (7/5 ^ 1), collapses and revivals do not 
occur, and the model exhibits more complex dynamical 
behavior [2^. Our interest here is in the recurrence prop- 
erties of the field observable (a+a) . For illustrative conve- 
nience, we take the atomic oscillator to be in the ground 
state and the initial state of the field to be either a CS 
or a PACS. We shall denote the corresponding states of 
the total system by |a;0) and |(a,m);0), respectively. 

In the case of weak nonlinearity, the return maps, 
recurrence plots, and first return distributions indicate 
quasiperiodicity of the field observables, remaining qual- 
itatively independent of whether the initial state of the 
field is a CS or a PACS, and of the location and size of 
the cell in the coarse-grained 'phase space' (i. e., the bin 
in (a+a)). These results are consistent with the occur- 
rence of revivals, and hence of returns of observables to 
their initial values. Figures [THTUl correspond to an initial 
state I a; 0), with v = 1 and 7/1? = 10~^. The return map 
in Fig. [7] and the recurrence plot in Fig. [5] reveal regular 
structures, while the distribution of the first return time 
for a typical cell is a sparse set of spikes. Figure [H depicts 
Fi for the cell C defined by the range [0.596,0.604] of 
(a+a). All these features are characteristic of quasiperi- 
odic dynamics. A numerical evaluation of the (invariant) 
density p of (^a^aj, the results of which are depicted in 
Fig. [TOl shows that this quantity is essentially made up 
of a large number of spikes. 

For strong nonlinearity (7/5 3> 1), however, revivals 
do not occur for generic conditions, and there arises a 
range of ergodicity properties, depending on the extent of 
coherence of the initial state of the field. For a significant 
departure from coherence (as measured by the ratio m/v 
for an initial PACS), or even for sufficiently large values 
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FIG. 8: Recurrence plot of (a^a) for the initial state |q;0) 
with u ^1- 7/g = 10"^ 
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FIG. 9: First return time distribution Fi of ^a^a) to the cell 
C, for the initial state \a; 0) with v — 1; j/g = 10~^. 
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FIG. 10: Invariant density of (a^a). Initial state \a;0) with 
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FIG. 11: Return map of the mean photon number (a^a) for 
the initial state |(a, 5); 0) with v — 5; j/g = 5. 



of v in the case of an initial CS, the dynamics of (a^a) ef- 
fectively becomes hyperbolic. The return map and recur- 
rence time statistics reveal this feature: the correspond- 
ing plots are in sharp contrast to those that arise in the 
quasiperiodic case. We illustrate this by means of plots 
for an initial state |(a, 5);0), with v = 5 and j/g — 5. 
The return map (Fig. Ilip and the recurrence plot (Fig. 
I12p no longer have well-defined patterns. The first return 
time distribution to a cell C = [12.455, 12.465], shown in 
Fig. I13[ approaches an exponential distribution (dotted 
lines), while the invariant density of (a^a), depicted in 
Fig. 1141 is consistent with a continuous distribution. 

As in the case of the single-mode model of Sec. the 
analysis of recurrence time distributions and recurrence 
plots is supplemented by a time series analysis leading 
to the estimation of the minimum embedding dimension 
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FIG. 12: Recurrence plot of (a^a) for the initial state 
|(q, 5); 0) with v = 5; 7/(7 — 5. 
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FIG. 13: First return distribution Fi of (a^a) for the initial 
state I (a, 5); 0) with u = 5; y/g — 5. 



P 



33 




3 




23 




2 




13 




L 




0^ 





11.S8 1.2.38 12M 



Increasing departure from coherence 
Initial state 





a ; 0> 


\(a,l) ; 0> 


(a, 5) ; 0> 




Y/g=IO^ 
|a|^v=l 


regular 


regular 


regular 


nlinearity 


Y/g=icr^ 

V=J 


Tf Crii laT" 


T"f CriilaT 


T"f CriilaT 


reasing no 


Y/g=5 
v=l 


regular 


regular 


chaotic 
exponential Fj 



c 




chaotic 


chaotic 


chaotic 


1 


Y/g=5 


exponential Fj 


exponential Fj 


exponential 




v=5 












chaotic 


chaotic 


chaotic 




Y/g=5 


exponential Fj 


exponential Fj 


exponential Fj 




v=10 






X ~ 1.0 

max 



TABLE I: Qualitative dynamical behavior of the mean pho- 
ton number of a single-mode electromagnetic field interacting 
with a nonlinear medium. 



and the calculation of the maximal Lyapunov exponent 
Amax- We have use d [2611 a robust algorithm developed 
by Rosenstein et alJA3^ and KantzM^ for the reliable 
estimation of Amax from data sets represented by time 
series. Table |T] gives a capsule summary of the repre- 
sentative results on the two-mode model of this Section. 
The characterizer 'regular' is used in those cases in vifhich 
Amax = to within numerical accuracy, while 'chaotic' in- 
dicates those corresponding to a positive value of Amax- 
It is immediately apparent that the latter is accompa- 
nied by an exponential recurrence-time distribution, in 
keeping with the known results for hyperbolic dynamical 
systems. As further corroboration we have also checked, 
in all these cases, that the distribution of two successive 
recurrences fits the next term in a Poisson distribution. 
The trends in the foregoing results are also in agreement 
with the results presented elsewherefisi on the entangle- 
ment properties of the model, including, in particular, 
the behavior of the subsystem entropy corresponding to 
the field mode. 

It is interesting that, even in as simple a model as the 
one considered here, where one of the components of a 
bipartite system effectively behaves like an open quantum 
system, the dynamical behavior can display such a range 
of ergodicity properties. The implications for genuine 
open quantum systems are therefore noteworthy. 



FIG. 14: Invariant density: Initial state 5);0) with v — 5; 
1/9 = 5. 
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